Spherical Collapse Model 
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We studied the spherical collapse model in the flat ACDM cosmology and provided exact and 
analytical formulaes for the calculation of the two important parameters in the application of Press- 
Schechter theory, the critical density contrast and the ratio of cluster/background densities at the 
virialization point in terms of a media variable which can be solved precisely by numerical methods. 
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INTRODUCTION 

Press-Schechter theory has been written into text- 
books, but two important parameters Sc and of it 
extracted from the studying of spherical collapse mod- 
els are not determined well enough. |^ are two early 
works studying the spherical models. They considered 
this model in a background universe without dark ener- 
gies, are four later works studying this model 
in the flat ACDM or open cosmologies. Our purpose here 
is to provide a simple and direct method for the calcu- 
lation of Scii^mQ, o-c) and Ad^mo, o,c)- Our background 
is flat ACDM cosmology but the method can be easily 
generalized to the QCDM as well as non-flat cosmolo- 
gies. Comparing with the existed results on this model, 
our results are analytical and exact. 



PRESS-SCHECHTER THEORY 

Press-Schechter Theory li| predicts that the fraction 
of volume which has collapsed at a certain red-shift z is 



f,ou{M{R),z) 



27ra{R, z) Js^ 
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Here, R is the radius over which the density field has been 
smoothed, <j{R, z) is the rms of the smoothed density 
field . <5c is the threshold of density contrast at time z 
(count stopping time) beyond which objects collapse. 

In existing literatures, there are two conventions for 
the definition of a{R,z) and Sc- In the first convention 
[ilSII3> '^(R^z) contains growth factor, while Sc only 
depends on the partition of the cosmological component 
today and the count stopping epoch ac- In the second 
convention j^, a{R) does not contain the growth factor, 
but Sc contains both the cosmological component parti- 
tion effects and the growth factor effects. We will take 
the first convention in this paper. Our Sc here will be 
corresponded to the Sc{0) of which was plotted in the 
upper panel of figure. 1 of it. 

Operationally, Sc is obtained by extrapolating the pri- 
mordial perturbation to the collapse epoch using the 



growth law of linear perturbation theory, i. e. 
,p,ncia,ac) 1 



Sc = 
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Where pmc and pmb are the mass density of the clus- 
ters and the background respectively; while Di (a) is the 
growth function of linear perturbation theory jj] , 



Di{a) 
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It can be shown that Di(a)a^a 
of M we can show that 



a. Using the method 
oc (1 - a • a) (4) 



So, 



Sc = 5a-Di{ac). (5) 

In a given cosmological model, Di{a) are known 0, the 
purpose of studying spherical collapse model is to deter- 
mine a. 

It should be noted that besides the partition of the 
cosmological components and the observational epoch, 
the value of Sc is also quite dependent on the choice of 
smoothing window used to obtain the dispersion a{R, z) 
0. We will not consider this effect in this paper. 

SPHERICAL COLLAPSE MODEL 

To calculate parameter a of eq®, let us start with the 
Friedman equations for both the over-dense region and 
back ground cosmology: 
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, where r and a denote the radius of the over-dense region 
and scale factor of background universe respectively; k 
is a constant; while subscript ta and c denote the turn 
around and collapse point time. Because at the turn 
aroimd time r = 0, 



SttG 



[Pmc,ta + PAo]- 
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Dividing eq© by (Uf and using © we get 



Pmc,tar ^rta + PAOr'^r^^ - {Pmc,ta + PAo) 



After letting 
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eqljnj becomes 
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We will not try to solve this equation, we derive it here 
because we need to use its asymptotic behavior to calcu- 
late 6c- 

In the ACDM cosmologies, the quantity C and ly in 
ea (|ll|) are only functions of ilmo and Oc- They do not 
depends on x or y. Using eql@}, we can write 



CHl-a-atax). (12) 
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So, 



dx 



[yx^o] = (1 - 2a • atax). 



(13) 



Substitute eas H12() and H13() into ea Hll|) . expand it and 
keep only the linear term, we get 



1 1 2 



(14) 



Substitute this result into eq©, we get: 
3 1 

5c(f^mO,ac) = -aiaV^^ + v ■ ataC^Di{ac) (15) 




FIG. 1: (5cS dependence on the cosmological models charac- 
terized by Q,mO and the collapse epoch a^. 




FIG. 2: dependence on Q.mO and the collapse epoch matter 
percentage Q.mb,c- Points on the same space curves have the 
same collapse epoch Sc- The outmost curve has Oc = 1, the 
next one has Oc = 0.9 and so on. The deeper a curve lies in 
the paper, the smaller Uc it has. 



This is the first important result of this paper. We 
will derive the differential equation for determining 
C,(^rnO,o,c) and the fitting formula ea l|29|l for it as well 
as the analytical solution for v{^rnOi ^c) ca (l32ll (and the 
notes there) in the technique details section. Comparing 
with the method of 0, our method here is more 
easier to operate and the result is more simply-looking. 

In FIG^ we plotted dependence on Q.raO and ac- 
From the figure we see that if Jlmo is fixed, then 5c is a 
decreasing function of ac- Physically this is because, in 
the given cosmological model, rimo fixed, the earlier an 
objects collapse, the more denser it is required to assure 
collapse occur successfully. On the other hand, if ac is 
fixed, then 5c increases as ilmo increases. Superficially 
this means that in those lower matter (thus high dark 
energy) percentage universes, collapse occurs more eas- 
ily! This is not the fact. Because in the lower matter 
percentage universes, density perturbation grows more 



slowly. As the result of linear extrapolation of the pri- 
mary perturbations eq|(2l, the more smaller growth fac- 
tor comes from the lower matter percentage suppresses 
the primary perturbations more strongly so gives more 
smaller 5c- 

To compare our results with that of Q , we plotted 5'^s 
dependence on Vlmo and Vlmb.c in FIGI21 From this figure 
we see that what was shown in the upper panel of figure. 1 
of 8] is the projection of our results in the 5c-^mb,c plane. 

To this point in this section, we imagined the evolu- 
tion of an over-dense region as: as long as it is over-dense 
enough, it will grow from an r = point, to the max- 
imum radius, then collapse to r = point. Factually, 
before the region collapse to the r = point, pressures 
from the random moving of the materials inside the re- 
gion will balance their self-gravity and the system will 
enter virialization status. At this point, the second im- 
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portant parameter Ac := in the application of 

Press-Schechter theory [EISIE^ appears. Assuming that 
at the coUapse point, the system has viriahzed fully, we 
can write: 



A _ Hmcyuta)' ta' c _ /- ' ta "c /-.f,-^ 
c ~ , \ -3 ~ s 3 ■ \^^) 

Pinb{ata)ataac '^c a-ta 

According to virial theorem and energy conservation law 
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we have 
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Looking as an equation for — , ea (|18|l can be solved an- 
alytically. Substituting the solution into ea H16|l . we get 



Ac(JlmO,ac) = C ■ 
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f(pt,P-c) = [{pt + 1) - -r^^^ 



+ {pt + l) + f{pt,Pc)]} (19) 
27 



= 7^ Pc = 



^ fic{27 fic-S{jit + m]^^^ (20) 
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Ea p9|) is the second important result of this paper. In 
ACDM cosmology, all the quantities on the right hand 
side of ea l|19|l are known functions of ^Imo and Uc- We 
will give the relevant formulaes ea (|28|l . (|29|l and H32(l in 
the technique section. Comparing with our results 
here is exact instead of approximated. 

To see the physical meaning of ea (jl9|l clearly, we plot- 
ted A^s dependence on fimo and Qc in FIG|31 From the 
figure, we see that Ac is an increasing function of Oc, i. 
e., the latter a region collapses, the more denser should it 
be. This can be understood physically. Because at more 
later times, dark energy's percentage will be more larger. 
To cancel its counter-collapse effects, more denser mat- 
ter is required to assure collapse. It is worth noting that 



FIG. 3: Left panel: A^s dependence on the collapse epoch ac, 
different curves denote different cosmological models indexed 
by flmo- Right panel: Aj.s dependence on the cosmological 
models, different curves denote different collapse epoch in- 
dexed by Uc- 
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FIG. 4: Right panel: A'^^^^f^p^s dependence on flmo and the 
collapse point Qrnb,c, the same curve has the same collapse 
epoch. Right panel: projection of three space curves from the 
right panel. 

used a different definition of Ac(ecf) := p^^^'^i (so 
different mass-temperature relations, please comparing 
eq(2.2) of Q and eq(4) of 0), we reproduce the results 
there in FIG0] From the right panel of this figure we see 
that, what was plotted in the lower panel of figure. 1 of 
is the projection of our results on the /S.c-^mb,c plane. 



TECHNIQUE DETAILS 

To calculate parameter ^ appearing in eas ljlSIl and 
()19|l . we write down Friedman equations for the back 
ground cosmology and the space-space component of Ein- 
stein equation for the radius of the over-dense region, 

-i:-{p^b + Ph)- (22) 
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(3pA -I- PA -I- Pmc)- 



(23) 



Using the notations introduced in ea (|lU|l . eas (|23|) and 
(I22|l can be rewritten as 

ST:Gpmb,ta 1 



y 



(24) 
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jwhere we preserved the equation of state coefficient w 
as a general parameter. It can even be a function of a or 
X. The above two equations can be translated into: 



dx 2 X 
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Ea H2()l) is a second ordinary differential equation, it con- 
tains a free parameter C,, but satisfies three boundary 
condition: 



— K^O — 
X 



Ja^O 



y\x=i = 1, y'\x=i = 0. 



(27) 



This is a two point boundary condition problem. It can 
also be looked as an eigen-value problem and solved nu- 
merically by the method described in §17.4 of 
Noting 



1 



:,(28) 



Substitute eqlj^U into H30() and solve it analytically, we 
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(32) 



where is today's ratio of dark-energy/mass density 
and V ~ i^oo-ta- m proposes that when eas (|^ and 
solved, using eq(A12) and (A13) of it to calculate 6c di- 
rectly. This method has difficulty to get high precision. 



CONCLUSIONS 

We studied the spherical collapse model in the flat 
ACDM cosmology and provided exact and analytical for- 
mulaes ea (|15f) and H19|l for the calculation of the two im- 
portant parameters Sd^mo, Oc) and Ac{flmO, o-c) in terms 
of the media variable C ea (|29|l . We reproduced the ex- 
isting results in the literature, but our method is more 
easier to operate and the result is more simply-looking. 
Our result will be useful in the studying of galaxy clusters 
evolution and the application of Press-Schechter theory. 



we see that C, solved from ea H2()|) only depends on ^mb,ta- 
So we can choose ata — 1 and solve ea (|26|l for different 
values of f2mO to get the appropriate C and fit the result 
as C(rimo)- Then change the function into C(^mfc,ta) to 
get C's dependence on and Oc. We do so by software 
Mathematica and confirmed the results by the method 
of and get 



.37r, 
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/-fa \ — N2o-0-7384+0.245ia„,t,| 



(29) 



Now the final question we need to answer is to express 
the scale factor of turn around point with that of the 
collapse point. Note that, in the unperturbed spherical 
cases, clusters' formation process is symmetrical about 
the turn aroimd time tta, so 



ic — 2tiQ. 

According to Friedmann equation 
the time-scale-factor relation is 



(30) 
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